Topological Qubits with Majorana Fermions in Trapped Ions 
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We propose a method of encoding a topologically-protected qubit using Majorana fermions in 
a trapped-ion chain. This qubit is protected against major sources of decoherence, while local 
operations and measurements can be easily realized. Furthermore, we show that an efficient quantum 
interface and memory for arbitrary multiqubit photonic states can be built, encoding them into a 
set of entangled Majorana-fermion (MF) qubits inside cavities. 
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The search for a decoherence free qubit is a major chal- 
lenge in the field of quantum computation [1 . Topologi- 
cally protected systems offer promising properties for the 
building of a fault-tolerant quantum memory [2 . How- 
ever, the realization of topological quantum memories up 
to now represents a challenging open problem. 

The discrete quantum wire model by A. Kitaev is one 
of the simplest systems supporting topological phases [3] . 
In this model, the signature for the topological nontriv- 
ial phase is the presence of unpaired Majorana fermions. 
Among their properties, these fermions coincide with 
their own antiparticles and support non-Abelian statis- 
tics Since its appearance, the model has attracted 
much interest, with diverse proposals for physical im- 
plementations, including superconducting heterostruc- 
tures f5^,^ and optical lattices [7j|8]. However, a clear ex- 
perimental evidence of unpaired Majorana fermions has 
not been achieved so far. Recently, it has been also sug- 
gested that superconductor-based Kitaev wires are prone 
to decoherence [9]. Moreover, building a qubit made out 
of Majorana fermions requires exceptional system con- 
trol. 

Trapped ions are highly controllable quantum sys- 
tems [10 . They can be cooled down to form crystals, 
easily initialized by optical pumping, manipulated with 
lasers, and efficiently measured. They offer one of the 
most reliable and scalable implementations for a quan- 
tum simulator [TTl [12]. Some examples are spin sys- 
tems [T3HI8], relativistic quantum mechanics [19-24 , 
quantum field theories [25], and fermionic systems [26 . 
A proposal for realizing another topological system, i.e., 
Kitaev honeycomb model, was put forward in Ref. [27] . 
However, the complexity of the honeycomb model, re- 
quiring three different kinds of interactions, XX, YY ^ 
and ZZ, is much higher than the one of the wire model, 
that only requires an XX interaction. Another proposal 
for topological systems in trapped ions [28], not involv- 
ing Majorana fermions, makes use of a highly nonlocal 
Hamiltonian. 

In this Letter, we propose the implementation of Ki- 
taev wire model [3] in trapped ions. By a mapping of 
Kitaev Hamiltonian onto a spin model, we show that 



this system can be realized in a trapped-ion chain with 
current technology, and a Majorana-fermion (MF) qubit 
can be encoded. This qubit will be topologically pro- 
tected against major sources of decoherence for longer 
times, constituting an efficient quantum memory. The 
proposed implementation, already valid for 3 ions, allows 
for the straightforward realization of local operations on 
the MF qubit and for an efficient readout of its state. 
Moreover, our proposal could allow for the first experi- 
mental creation and test of Majorana fermions, and the 
study of this exotic phase of matter. We also show that 
a quantum interface between highly entangled incoming 
photonic states and MF qubits can be implemented by 
grouping many of these basic units. To this end, we sug- 
gest the use of an array of ion chains inside a set of cav- 
ities as a possible experimental realization. 

We begin by considering the Kitaev Hamiltonian [3] 
made up of N spinless fermionic sites, 

N-l 

-/i^ (a]aj - - j , (1) 

where w is the hopping energy, /i is the chemical potential 
and A is the order parameter of the pairing term, which 
mimics a p-wave spinless superconductivity. Given that 
A = I Ale*"^, one can define a set of 2N Majorana fermions 
C2j-i = e^2aj +e *2aj^ C2j = —ie'"^aj + ^2aj^ and 
the Hamiltonian becomes 

Hm = [-/iC2j-lC2j + + |A|)c2jC2j+l 

j 

+ (-«;+|A|)c2,-iC2,-+2]. (2) 

Under the parameter regime ja = 0^ \ A\ = w > the Ma- 
jorana fermions ci and C2n disappear, and Hm becomes 
diagonal in the basis aj = ^(c2j +ic2j+i), j = 1,...,X — 1. 
One can pair the two outer Majorana fermions into an 
additional complex fermion aj^ = (ci — zc2Ar)/2. The 
ground state of the Hamiltonian under this regime is 
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spanned by the occupancy subspace with respect to this 
nonlocal fermion, and is separated by a gap 2w from the 
excitation of the system. The ground state is two-fold de- 
generate, and is spanned by the states |^o) and de- 
fined by ai|^o,i) = 0, Vz = l,2,...,Ar-l and aul'^o) = 0, 
l^i) = <^mI^o)- It was proved [3 that these states sur- 
vive in the regime < 2w. Under this regime one has 
a nontrivial value for the topological invariant Z2, that 
labels the two different topological phases of the system. 

Trapped ion implementation.- In order to implement 
Kitaev model in an ion string, we take into account that 
Hamiltonian ([T]), for the parameter regime (/> = 0, which 
does not reduce the generality of the model, and |A| = 
w > 0, can be mapped onto the transverse field Ising 
model by using a Jordan- Wigner transformation [29l [30] , 

N-l N 

where J = w is the exchange coupling and hz = — /i/2 
is a transverse magnetic field. The mapping of the 
Majorana fermions into the spins is C2j = Ilti^fe^j^ 

C2j-i = ni=i ^k^j' definition of aM, ^i, and |^o,i) 
previously defined change accordingly. The two ground 
states 1^0,1) of Hg in Eq. ([3| consequently span a sub- 
space that is protected from higher-level excitations by 
the energy gap 2w = 2 J. 

Hamiltonian Q can be implemented with standard 
trapped- ion technology [iO\. We consider a ID string 
of two-level ions coupled through the joint motional 
modes by means of external lasers. One possibility is to 
apply Raman lasers for the spin-spin interaction and for 
the magnetic field at the same time [13l [Mj [3Tj . By lo- 
cally addressing each ion with different Raman beatnotes, 
in a multiple M0lmer-S0rensen configuration, one can 
achieve the Ylf=i ^j^±4-i interaction with open bound- 
ary conditions in Eq. ([3|, while a single laser can imple- 
ment the a J part ^31^ . An alternative possibility is to 
apply a Trotter decomposition, alternating between the 
different interactions pT} [TS]. 

Here we are interested in the pure Majorana regime 
in which J ^ \hz\. The aim is to achieve this regime 
in the always-on interaction, such that the degenerate 
ground state of Hamiltonian ([3| encodes the topologi- 
cally protected subspace. This degenerate ground state 
is composed of two Greenberger-Horne-Zeilinger (GHZ)- 
like states [32 , that are highly entangled, e.g., for odd 
TV, 1^0) = (I ... H - I ... and 

1^1) = (I ... ^) + I ... ^))/V2, while for 

even TV the relative signs of the two terms in both |^o) 
and 1^1) are interchanged. These states have different 
parity P = - JJi^i^ ^l^o,i) = ±|^o,i), which is a con- 
served quantity with respect to Hamiltonian (|3|. To ob- 
tain them, one may use standard quantum computation 
techniques [32 . Indeed, GHZ states of up to 14 qubits 
have been obtained experimentally with trapped ions. 



On the other hand, they can also be easily achieved us- 
ing an adiabatic evolution [33 from the opposite regime 
in which J <C \hz\. The ground state of this Hamil- 
tonian -hzT^j^j for hz < 0, l^lf^*) = ... ;) is 
mapped adiabatically onto the ground state |^o) of the 
Hamiltonian —J^jSi cr^cr^j^^- In turn, one specific lin- 
ear combination of the first excited states of —hz 
1^1"^*) = Xli \ ... ti ... \)i is mapped onto the ground 
state 1^1) of —JYl^^i I^or instance, for = 3, 

Ci = C3 = 1/2, and C2 = 1/a/2. Notice, in addition, 
that the two states will not be mixed during the adia- 
batic evolution due to their different parity. In order to 
create state one may just apply standard optical 

pumping. In order to create |^\"^*) one may apply an in- 
homogeneous Tavis-Cummings Hamiltonian [34] with ap- 
propriate couplings, engineered with an inhomogeneous 
red-sideband interaction upon the ions, using the center- 
of-mass mode (e.g., with displaced equilibrium positions 
in the trap by a tailored potential [35]). 

Quantum interface.- This adiabatic mapping may be 
most appropriate for transferring an arbitrary qubit state 
that comes into the system, e.g., a photonic qubit, to 
the MF qubit. One can begin with a photonic incom- 
ing state, a|0) + /^|1), where |0), |1) are Fock states: by 
using a quantum interface with a cavity, state |0) can 
be mapped to and |1) to by a collective 

excitation [36 . We point out that in order to get high fi- 
delities in the coupling to the incoming photon, chains of 
around 20 ions at least should be employed [36 , though 
for the MF qubit itself, as we will show, 3 ions are in 
principle enough. Subsequently, the adiabatic transfer 
will produce the same qubit superposition a|^o) 
in the MF qubit. An extension could be considered with 
K copies of this system, with a highly-entangled incom- 
ing photonic state, ^ Q^...^^ Here each photonic 
qubit is mapped similarly to the corresponding MF qubit, 
giving i.e., an arbitrarily entangled 
state of nonlocal Majorana-fermion qubits. This could 
be used as an efficient quantum memory. In Fig. ^ we 
show a scheme of the proposed setup. 

Local operations.- One can efficiently implement lo- 
cal operations upon the MF qubits with our setup. 
In general, these local operations upon these nonlocal 
MF qubits translate into nonlocal operations upon the 
physical qubits, i.e., the trapped ions. On the other 
hand, these nonlocal spin operations upon the N ions 
of the qubit can be efficiently implemented, as recently 
shown [W, '37] , in terms of a reduced number of global 
and local laser interactions. More generally, the complete 
set of local operations, i.e. Pauli matrices, upon the non- 
local MF qubit, cr^pQ, cr^pQ' ^mfQ' t)e related 
to the single-ion Pauli operators {cr^'^'^}i=i,..,Ar through 
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FIG. 1. (a) Array of ID ion crystals inside K optical cav- 
ities for the quantum interface between a multiqubit en- 
tangled photonic state and an MF qubit multipartite state. 
The shaded region denotes a single MF qubit. (b) Fidelity 
F = |(^o|^(t))|^ for the evolved state |V^(t)) from the initial 
state r^o^^*) using the adiabatic transfer protocol with Hs in 
Eq. (|3]), for AT = 3. 
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FIG. 2. (a) Energy splitting F between |^o) and as a 

function of number of sites N ^ in units of J = 1, computed nu- 
merically for hz — 10~^J and random nearest-neighbor cou- 
plings Ji,i+i with \ J—Ji^i^i \ < O.OIJ. (b) Survival probability 
F = |(^o|V^(t))|^ for the evolved state |V^(t)) from |^o) un- 
der the dynamics of Hamiltonian (|3| (with hz — 0) (solid) or 
without it (dashed) plus random local operators of modulus 
10""^ J in direction Z, proportional to erf, for N — ?>. 



the expressions 

^MFQ = ajvf + ttM = / (8) / (8) ... / , 

^MFQ -^(^M - ^m) = -CTn ^ 8) ... , (4) 

^MFQ «M^M - ttM^M = ^ ^ 8) crf (g) cr^. 

As can be appreciated, cr^pQ are highly nonlocal oper- 
ations upon the ions. Nervertheless, they can be im- 
plemented efficiently with a reduced number of lasers as 
recently shown [26l[37]. Thus, with our setup we have a 
fully controllable Majorana-fermion qubit. 

Measurement.- Another main appeal of our proposal 
for encoding MF qubits in trapped ions is the ease of per- 
forming measurements on this system. A projective mea- 
surement upon the basis {|^o), which is equivalent 
to detecting the parity operator P = — Fli^f^ amounts 
to one local measurement of operator per ion. This 
is the easiest detection performed in trapped ions and 
can be done, with standard resonance fluorescence, with 
fidelities larger than 0.99 [10 . Combined with the lo- 
cal operations exposed above, this allows one for the full 
tomographic reconstruction of the MF qubit. 

Error and decoherence protection.- This implementa- 
tion, as has been shown above, contains a degenerate 
ground state, {l^o)^ l^i)} which is protected by a gap 
from higher-level excitations. Observation of this gapped 
topological system will be in itself interesting. A conse- 
quence of this is that local operations af'^, which couple 
ll^o), l^i)} with higher energy states, are topologically 
supressed. Thus, magnetic field fluctuations along these 
directions will not have an effect upon the system. On 
the other hand, local af operations realize swap gates 
between the |^o) and states. In our trapped ion 
setup, random ambient magnetic fields along X direction 
will rotate fast in the interaction picture with respect to 
the trapped-ion qubit energy in which Eq. (|3| is com- 
puted [31], such that their contribution to decoherence 
will be negligible. In general this will also happen with 



the F-direction ambient magnetic fields, such that these 
will be doubly protected: by the topological Hamiltonian 
and by the single-ion qubit energy. 

Feasibility of the implementation. - The qubit encoded 
in Kitaev wire model [3 is robust in the large TV limit, 
for the parameter range < 2w. On the other hand, 
for finite the degree of imperfection in the proto- 
col depends on how much the system deviates from 
the perfect parameter regime w = |A|, /i = 0. In- 
deed, the appearing energy splitting F between |^o) 
and l^i), which breaks the ground-state degeneracy and 
qubit coherence, is of the order of F oc exp(— A^/no), 
where n^^ = min{| In | In and x± = (— /i ± 

-4w^ ^4\A\^)/[2{w + |A|)]. We consider J - 6 
kHz for the 3 ion case. For realistic imperfections in the 
Rabi frequencies of the lasers of 1%, that induce the same 
order of imperfection in = A = J, and magnetic field 
fluctuations in /i = —2hz ~ 2x 10~^ J, we have no — 0.14. 
1% site-dependent deviations | J — J^^^+i|/J will slightly 
increase this value (see below), but even for reduced 
F will be negligible. This makes the implementation of 
the wire model and Majorana fermions in trapped ions 
a real possibility with current technology. We show be- 
low that considering a linear chain of just A^ = 3 ions, 
one may encode a topologically-protected qubit subspace 
with very low decoherence. 

In order to test the feasibility of the implementation 
with realistic trapped-ion systems, we realized numerical 
simulations with different kinds of imperfections. With 
respect to the adiabatic protocol with Hs in Eq. (|3|, 
we have computed the fidelity F = K^olV^l^))!^ for the 
evolved state |'^(t)) from |^q"^*), making the evolution 
\hz\/J > 1 ^ \hz\/J < 1, for A" = 3, see Fig. [i]d. We 
considered the ideal case with no spurious hz final value 
and equal nearest-neighbor couplings Ji,i+i = J, and the 
diabatic error is of 10~^. Introducing realistic random 1% 
deviations of Ji,i+i from J due to laser-intensity fluctua- 
tions and 10~^ J final fluctuations in hz due to local ambi- 
ent magnetic fields did not modify appreciably the proto- 
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col. We plot in Fig. [2] (a) the energy splitting F between 
1^0) and 1^1) as a function of A^, in units of J = 1, com- 
puted numerically for hz = 10~^ J and random nearest- 
neighbor couplings Ji,i+i with \ J — Ji,i+i| < O.OIJ. The 
scaling is exponential, as expected. For A/" = 3, we have 
T 1.6 X 10~^ J, i.e., a negligible dephasing error. Thus, 
this topological encoding will in principle work already 
for 3 ions, which is feasible with current technology. 

To show the subspace topological protection of the 
degenerate ground state {l^o), l^i)} with respect to 
coupling to outer states, we plot in Fig. [2] (b) the 
survival probability K^olV^l^))!^ for the evolved state 
\ip{t)) from 1^0) under the dynamics of Hamiltonian (|3| 
(solid) or without it (dashed) with random fluctuations 
Shi ^ ^^~^J of for = 3. This is the kind of local 
operators that couple spanjl^o), l^i)} with outer states 
to the subspace. It is clearly appreciated that the survival 
probability inside the subspace is significantly increased 
upon evolution with (|3| . To estimate the order of magni- 
tude of the protection, we consider slowly- varying, local, 
random magnetic fields along af . We remark that this is 
the main fundamental decoherence source in most exper- 
imental trapped ion setups [38 , with coherence lifetimes 
of about 30 ms. We have that states l^o)^ will couple 
to each other or to outer states only through higher-order 
processes in perturbation theory. The largest contribu- 
tion to decoherence is at most of order ~ Sh^/Ag. Here 
Shz is the average local magnetic field perturbation, that 
we take, as in the numerical simulation, to be 10~^J, 
and Ag = 2J is the gap between the topological ground 
states 1^0)7 1^1) and the excited states. For this we have 
assumed that these spurious magnetic fields change in 
time much more slowly than the frequency of the gap. 
Accordingly, the effective Rabi frequency of the error is 
about 5 X 10~^Shz^ i.e., several orders of magnitude re- 
duced with respect to the case without topological pro- 
tection [without evolution by ([3|], that is of order Shz. 
Considering that motional heating rates in trapped ions 
can be of a few phonons per second and spontaneous 
emission lifetimes of more than one second [38 , with this 
proposal one may improve coherence lifetimes by more 
than one order of magnitude with current technology. 

Conclusions. - We showed that trapped-ion chains can 
implement a topologically protected qubit subspace by 
means of Majorana fermions. We proved that MF qubit s, 
starting already with 3 ions, can outperform the usual 
ionic qubit coherence time by more than one order of 
magnitude, yielding an efficient quantum memory. More- 
over, a quantum interface with photonic states could be 
implemented with longer chains, allowing us the realiza- 
tion of two-qubit gates among the MF qubits and quan- 
tum communication. Our proposal may also give rise to 
the unprecedented creation of Majorana fermions in a 
controllable quantum platform. 
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